Coherent Josephson qubit suitable for scalable quantum integrated circuits 
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We demonstrate a planar, tunable superconducting qubit with energy relaxation times up to 
44 /is. This is achieved by using a geometry designed to both minimize radiative loss and reduce 
coupling to materials-related defects. At these levels of coherence, we find a fine structure in the 
qubit energy lifetime as a function of frequency, indicating the presence of a sparse population of 
incoherent, weakly coupled two- level defects. This is supported by a model analysis as well as exper- 
imental variations in the geometry. Our 'Xmon' qubit combines facile fabrication, straightforward 
connectivity, fast control, and long coherence, opening a viable route to constructing a chip-based 
quantum computer. 



One of the outstanding challenges in building a quan- 
tum computer is to balance coherence, connectivity and 
control in the qubits. Superconductivity provides an ap- 
pealing platform because it allows for scalability: the 
conduction electrons condense into a macroscopic quan- 
tum state, and large quantum integrated circuits can be 
made with many elements having individual control lines. 
However, quantum coherence in superconducting circuits 
has proven to be very delicate, as it is easily disturbed 
by material defects, electron system excitations, and ra- 
diative coupling to external wiring pj-l2|]. To minimize 
these and other effects, many groups have recently begun 
embedding qubits in three-dimensional superconducting 
cavities. These 3D qubits show high coherence, with en- 
ergy relaxation times in 3D transmon qubits between 30 
and 140 /is %%. 



Here, we demonstrate a new design for a fully pla- 
nar superconducting qubit, based on the planar trans- 
mon [lO|, UjJ, with energy coherence times in excess of 
40 /is. Our approach balances coherence, connectivity, 
as well as fast control. We design the qubit with high- 
quality coplanar waveguide capacitors, motivated by the 
recent advances with superconducting resonators, yield- 
ing a modular design with straightforward connectivity. 
The qubits are frequency- tunable, which allows the im- 
plementation of fast two-qubit gates: a controlled- Z gate 
[12h14| can then be implemented with high fidelity in 
25 ns [l5|. With the coherence time exceeding single and 
two-qubit gates times by three orders of magnitude, we 
believe that our device provides a key ingredient for im- 
plementing a surface code quantum computer [16]. 

We also find that the small remnant decoherence in 
these qubits comes from a sparse bath of weakly cou- 
pled, incoherent defects. These defects are clearly visi- 
ble in the measured quantum time-resolved spectroscopy, 
and give rise to frequency-dependent variations in the en- 
ergy relaxation rate. These results may also explain the 
variations observed in lifetimes of 3D transmon qubits. 
We elucidate this physics and improve the coherence by 
varying the capacitor geometry. 



Our device is shown in Fig. [T^, formed by pattern- 
ing the Al metal (Hght), exposing the sapphire substrate 
(dark). The qubit is the cross-shaped device. Its four 
arms connect to separate elements, each having a differ- 
ent function: a coplanar waveguide resonator for readout 
on the top, a quantum bus resonator on the right to me- 
diate coupling to other qubits, XY control on the left to 
excite the qubit state, and Z control on the bottom to 
tune the qubit frequency. The cross is the qubit capaci- 
tor, which connects at the bottom to the tunable Joseph- 
son junction, formed by the rectangular ring-shaped su- 
perconducting quantum interference device (SQUID), see 
Fig. [TJd. The rectangular ring is intersected by two iden- 
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FIG. 1: (Color online) (a) Optical micrograph of the planar 
'Xmon' qubit, formed by the Al superconducting film (light) 
and the exposed sapphire substrate (dark). The qubit is ca- 
pacitively coupled to a quarter wave readout resonator (top), 
a quantum bus resonator (right), a XY control line (left), and 
inductively coupled to a Z control line (bottom). The Xmon 
arm length is L. (b) The inset shows the shadow evaporated 
Al junction layer in false color (blue). The junction size is 
0.30 X 0.20 /xm^. The capacitor central linewidth is 5*, the 
gap width is W . (c) The electrical circuit of the qubit. 



tically sized Al tunnel junctions (blue in Fig. ^). The 
electrical circuit is equivalent to that of a grounded trans- 
mon [lO| , with the capacitor in parallel with the tunable 
junction (Fig. [it). In a clear departure from the tradi- 
tional floating transmon with an interdigitated capacitor 
(IDC) [ul, we chose to form the qubit capacitor by in- 
tersecting two coplanar waveguide lines. 

In prior work, we showed that highly coherent coplanar 
waveguide resonators can be fabricated, having quality 
factors of about 1.5 x 10^ at the single photon occupa- 
tion level. These resonators were made from molecular 
beam epitaxy (MBE) Al on oxygen-cleaned sapphire [17| . 
This shows that a straightforward path to high coherence 
comes from a combination of I) MBE Al as high quality 
material, II) coplanar waveguides having low radiative 
loss, and III) embedding in a groundplane. We therefore 
embed the qubit in an uninterrupted groundplane, with 
thin Al lines at the capacitor ends tying the groundplanes 
together; this suppresses parasitic slotline modes in the 
control lines and resonators as well. 

Connectivity is accomplished by coupling each of the 
qubit 's arms to a distinct element with specific function- 
ality. Three of the connections are easily made with a 
coupling capacitor, as the qubit is connected to ground. 
An advantage of this approach is that each coupling can 
be individually tuned and optimized. To this end, we 
have also separated out qubit control. The XY control 
drive line is connected with a coupling Cc = 60 aF, which 
allows us to excite the qubit state in 10 ns but hardly af- 
fects coherence, with an estimated Ti of 0.3 ms. The Z 
control also combines speed and coherence. The drive 
line is galvanically connected to the SQUID to allow for 
a large inductive coupling with a mutual inductance of 
M = 2.2 pH. We are able to rapidly detune the qubit on 
the timescale of a nanosecond [18]. The measured para- 
sitic coupling between the Z line and the qubit gives an 
estimated Ti of ~30 ms [19]. 

The qubit design outlined above balances coherence, 
control and connectivity in the qubit, using a robust de- 
sign for the capacitor and a modular separation of func- 
tionality. Our approach differs on key points from previ- 
ous implementations. With this experimental nature in 
mind, we name our qubit the 'Xmon'. While the cross- 
shaped qubit capacitor may emphasize this name, more 
arms can be added to allow for more connectivity. 

We find a dramatic increase in Xmon energy coherence 
compared to the traditional planar transmon, measuring 
decay times up to Ti = 44 /is, see Fig. [2^. The qubit 
Ti is measured by exciting it with a 7r-pulse, and mea- 
suring its state after a variable time Ar. We find that 
the excited state probability decays exponentially. We 
find Ramsey and spin echo phase coherence times up to 
T2* = 15 /is and T2 = 20 /is at the flux insensitive point, 
respectively, see Fig. [JJd. The energy decay at this point 
is measured to be Ti = 18 /is. The dephasing envelopes 
are measured using tomography. The first pulse is X^^2^ 
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FIG. 2: (Color online) (a) Qubit energy decay at three nearby 
frequencies. The qubit frequency is adjusted by applying a 
rectangular pulse with length Ar on the Z line. The pulse 
sequence is shown in the inset, (b) Ramsey T2* and spin echo 
T2 dephasing envelopes at the flux insensitive point, measured 
by phase tomography. The inset shows the pulse sequence; for 
the spin echo we apply a refocussing pulse (dashed). We apply 
four phases to the last pulse for phase tomography, to measure 
the decay envelope. Spin echo measurements are limited by 
electronics to 14 /is. The energy coherence for the dephasing 
measurement was Ti — 18 /is. 



followed by X^/2^ -^-7r/2 7 ^7r/2 7 or ^-7r/2 (sGC iusct for the 
pulse sequence), producing fringes with different phases. 
The dephasing envelopes follow an exponential decay. As 
the limit of T2 = 2Ti has not been reached [10], this in- 
dicates the presence of an additional dephasing channel. 
This channel, as well as dephasing away from the flux 
insensitive point, is presently under investigation. 

The qubits used here had ground to excited state tran- 
sition frequencies around 6 GHz when unbiased, nonlin- 
earities around 230 MHz, and a ratio of Josephson to 
charging energy Ej/Ec ^ 95. We employ a dispersive, 
high-power single-shot readout scheme with a 70-85% fi- 
delity [20]. The readout resonator frequencies used are 
6.4-6.7 GHz, the loaded quality factor is Qi = 10^, and 
the resonator-qubit coupling strength is approximately 
40 MHz. Measurements were done in a dilution refriger- 
ator with a base temperature of 30 mK, with multistage 
infrared shielding [21]. Magnetic fields were reduced by 
room temperature and cryogenic magnetic shields, with 
nonmagnetic microwave connectors [22]. 

The results in Fig. [2] show that tunable superconduct- 
ing qubits with a planar geometry can have Ti values in 
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FIG. 3: (Color online) (a) Qubit spectroscopy for device with 
*S',VK=8,8 /im. A smooth curve is formed by the high trans- 
mission (grey), measured on resonance with the readout res- 
onator, which indicates when the qubit is excited, (b) Swap 
spectroscopy of the same qubit. The qubit is detuned from 
4 to 6 GHz (stepsize 2 MHz), and the delay time is varied 
from 100 ns to 150 /is. See the inset of Fig. [2^ for the pulse 
sequence, (c) Qubit relaxation rate at qubit frequencies from 
5.25 to 5.75 GHz, extracted from the data in (b). Sharp 
peaks above a background are observed. The peaks are fitted 
to Eq.[T] (solid lines). 



excess of 40 /is. In fact, this Ti corresponds to the MBE 
Al resonator quality factors [17|, for which Ti = Q/uj is 
also about 40 /is. The combination of long energy and 
phase coherence times compare well with previously re- 
ported values for planar superconducting Al qubit s: for 
transmons Ti = 9.7 jis and T2* = 10 us [23|, for charge 
qubits Ti = 200 jis and T^ = 0.07 //s [2i|, for flux qubits 
Ti = 12 /is and T2* = 2.5 /is |25], and for the fluxonium 
Ti = 10 /is and T2* = 2 /is [26]. In fact, the Xmon ap- 
proaches the long coherence found in 3D transmons [8|, |9| . 
Very rece ntly , TiN planar devices have shown long coher- 
ence [27, 28], encouraging using Xmon geometries with 
this material. 

We find that the energy relaxation depends on qubit 
frequency, as we observe different exponential decay rates 
for small changes in frequency. As shown in Fig. [2^, we 
find Ti values from 34 to 44 /is in a 4 MHz band near 
5.4 GHz. In order to elucidate this further, we performed 
a spectroscopic scan on the qubit, shown in Fig. [3^1. 
The qubit frequency displays the expected dependence 
on applied flux <l> [10|, varying smoothly without visible 
splittings, indicating that strongly coupled defects, which 
manifest as avoided level crossings [1], are virtually ab- 
sent. We then performed a quantum analogue of time- 
resolved spectroscopy (swap spectroscopy [29]), shown in 



50 

40 

30 
-w 20 



S=8,W=4 







6 4 




5 
f(GHz) 



FIG. 4: Frequency dependence of Ti for six qubits with differ- 
ent S and W (see Table [J). The frequency stepsize is 5 MHz 
for S,W — 8,4 /xm and 2 MHz otherwise. See Supplementary 
Material [ly] for the corresponding decay rates. 



Fig. [3b. The probability of the excited state (color) is 
plotted for Ar from 100 ns to 150 /is (logarithmic ver- 
tical scale) and qubit frequencies from 4 to 6 GHz. We 
find that the probability decays exponentially, but with 
a fine structure of variable energy relaxation, and dis- 
tinct peaks in the energy decay rate (Fig.[3fc). We do not 
observe any chevron interference patterns [29], implying 
no defects interact coherently with the qubit state. Af- 
ter cycling the temperature to 4.2 K the fine structure 
is altered, but the overall image remains unchanged. We 
count approximately 30 regions with reduced coherence 
(Ti < 8 lis) per GHz in Fig. [3b. 

We explored the dependence of the qubit coherence 
time on capacitor geometry, using six different qubit ca- 
pacitor designs; the central line width S^ gap width W 
and arm length L were varied, while the capacitance 
value [30| and junction parameters are kept the same. 
The parameters are listed in Table [T] see Supplementary 
Material [18] for a micrograph. We find that the swap 
spectroscopy measurements of the different designs share 
the same characteristics as shown in Fig. [3b: a fine struc- 
ture with varying exponential decay. The energy relax- 
ation times extracted from the measurements are shown 



TABLE I: Geometric parameters for the Xmon qubit capaci- 
tors as defined in Fig. [T] along with their frequencies. Groups 
of three qubits indicate that the devices are on the same chips. 



S (/im) 


8 16 24 


8 16 24 


W (/im) 


4 8 12 


8 16 24 


L (/im) 


130 130 130 


165 165 165 


/lO,max (GHz) 


6.094 6.158 6.071 


6.080 5.883 5.846 


nonlinearity (MHz) 


224 228 222 


220 225 223 



in Fig.|4j With larger capacitor size the overah energy re- 
laxation time increases: when changing S^W from 8,4 /am. 
to 16,8 /xm and 24,12 /im, the Ti improves from a band of 
values between 8-15 to 10-20, and 20-40 /is, respectively. 
Importantly, both the upper as well the lower bounds 
on Ti increase with capacitor width. This is repeated in 
the qubits with S^W ranging from 8,8 /im to 16,16 /im 
and 24,24 /im. The reduction of Ti at frequencies ap- 
proaching 6 GHz is due to Purcell decay into the readout 
resonator [Sjj. We emphasize that these Ti values are 
obtained in multiqubit chips including complete sets of 
control wiring. 

The improvement of Ti with increasing width is con- 
sistent with previous experiments on superconducting 



resonators [32|, |33|. Loss arises from the electric fields 
coupling to two- level systems with dipole moments [34], 
which reside predominantly in surface oxides and inter- 
faces. This loss depends on the participation ratio, which 
depends on the electric field distribution [35| . Widening 
the capacitor reduces the participation of the surfaces 
and thus the loss, a natural explanation for the approxi- 
mately linear increase in average Ti with width in Fig. HI 
On the other hand, the peaks in the decayrate are remi- 
niscent of experiments with phase qubits [1| , where local- 
ized features in the frequency dependence occur when the 
qubit couples strongly to two- level defects, often giving 
rise to splitting of the qubit frequency and the chevron- 
shaped signature of coherent swapping. However, the 
exponential decay in the Xmon qubit, with no signatures 
of swapping or splitting, suggests a different energy re- 
laxation mechanism. 

Here, we show how surface defects near the metal edges 
of the capacitor provide a natural explanation for the 
peaks in the energy decay. The key point is that loss 
arises from the qubit interacting with a sparse bath of 
incoherent defects, as indicated by our data: the sharp 
frequency dependence as well as the changes in fine struc- 
ture when cycling to 4.2 K are both clear indicators of de- 
fects. The absence of chevrons and qubit frequency split- 
tings correspond to incoherent interaction. The lower 
and upper bounds of Ti increasing with capacitor dimen- 
sion indicates the defects reside in the capacitor. 

We model a quantum system consisting of a qubit, with 
a frequency-independent background loss rate Fi^q and 
pure dephasing rate F^^q, and a single two- level defect 
with decoherence rate Fi^d and dephasing rate F^^d (see 
Supplementary Material |l8|). When the defect deco- 
herence rate exceeds the coupling strength ^, coherent 
swapping vanishes and an incoherent, exponential decay 
appears. From a two-spin Hamiltonian (see Supplemen- 
tary Material) , we derive the qubit energy relaxation rate 
Fi (in the limit Fi^d > 9 > Ti^q) to be 
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(1) 



Hence, each uncorrelated defect adds a single Lorentzian 
to the energy decay rate. We can roughly estimate the 
coupling strength g for a surface defect with dipole mo- 
ment p ~ 1 D at a distance x away from the metal edge. 
With the electric field given by E = B/^/x [36), and 
B from numerical simulations [37|, we arrive at g/27T ^ 
0.1 MHz {g = pE) for a defect located at x = 3 nm away 
from the metal. We apply our model to Fig. [3h and find 
that the peaks in decay rate can be described by a set 
of Lorentzians, with 1/(Fi,d/2 + F^^d) ~ 50 — 100 ns, 
consistent with defect decay rates measured in similar 
systems [ll, |38|, and with g/27r > 0.2 MHz, agreeing with 
incoherent loss. 

We can also estimate the number of individually resolv- 
able defects using two-level system physics developed for 
junctions. The substrate- metal interface in our devices 
was thoroughly cleaned [17], hence we assume that the 
bulk of strongly-coupled defects resides in the metal- and 
substrate-air interfaces, as they have the highest partici- 
pation ratios ^] . The defect density for AlO^^ in tunnel 
barriers has been established in measurements with phase 
qubits [l| , with the distribution over dipole moment given 
by po\/l ~ P'^/Pmax/P^ with po ^ 10^//im^/GHz, and the 
maximum dipole moment Pmax = 6 D. We take these 
numbers as representative and assume a 3 nm thick di- 
electric layer with defects [39]. The number of defects 
with coupling strength greater than gmin is then given by 
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with detuning A, and F = Fi,d/2 + F0,d + Fi,q/2 + F0,q. 



with the unit step function and E{r) the electric 
field at position r. Simulations using Eq. [2] as well as 
Monte Carlo simulations indicate N ^ 30 — 50/GHz, for 
gmin/^TT ~ 0.2 MHz. We emphasize that the simula- 
tions connect ^min to A^ with values which are close to 
what is observed experimentally. The Monte Carlo sim- 
ulations (see Supplementary Material [18]) also indicate 
that the bulk of strongly coupled defects reside within a 
~100 nm distance from the etched metal edges, where 
the electric fields are largest. In addition, we have simu- 
lated the qubit decay rate (see Supplementary Material), 
reproducing the experimentally observed features: the 
simulated decay rate shows both the peaks as well as the 
variation in the background. 

The good quantitative comparison between model and 
experiment gives compelHng evidence that a sparse bath 
of incoherent defects plays a major role in loss in highly 
coherent qubits. Our results may also explain previously 
reported anomalous behavior in planar transmon qubits 
with long coherence, for which the Ti has been reported 
to vary significantly between qubits, even on the same 
chip [23, 40]. This is consistent with a sparse bath of 
incoherent defects limiting the coherence, as in Fig. [3l 

In conclusion, we demonstrate energy coherence times 
exceeding 40 /is in tunable, planar superconducting 
qubits. We have achieved this using a geometry with 



low radiative dissipation and high quahty materials. At 
these high coherence levels we identify a novel decoher- 
ence mechanism, loss from a sparse bath of incoherent 
defects, which is apparent in the swap spectroscopy. We 
show how enlarging the capacitor improves coherence. 
Our qubits combine long coherence, easy inter connectiv- 
ity, and fast control, providing a key ingredient for the 
implementation of an on-chip surface code quantum com- 
puter. 
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DEVICE FABRICATION 

The devices were made in a two-step deposition pro- 
cess. The qubit capacitor, groundplane, readout res- 
onator and control and readout wiring were made in a 
first, separate deposition step. We used molecular beam 
epitaxy (MBE) Al deposited on a c-plane sapphire sub- 
strate. The Al film thickness is approximately 100 nm, 
deposited at room temperature. The sapphire substrate 
was cleaned by load- lock outgassing at 200 °C, followed 
by heating to 850 °C in ~ 10~^ Torr activated oxygen, 
identical to the process outlined in Ref. [l|. The first Al 
layer was patterned by a BCI3/CI2 reactive ion etch. 

In the final step, 0.30 x 0.20 fim^ Al tunnel barriers 
(30 nm bottom and 100 nm top layer thickness) were 
made using double-angle shadow evaporation. We used 



a high vacuum electron-beam evaporator, with a base 
pressure of approximately 5 x 10~^ Torr. We used the 
Dolan bridge technique with a poly (methyl methacry- 
late) /copolymer resist bilayer (approximate thickness: 
0.30 and 0.50 /im, respectively), patterned with electron 
beam lithography. In order to make galvanic contact 
between the first Al layer and the junction layer, we 
used a 3 min long Ar ion mill (beam: 400 V, 21 mA; 
beam width: ~3.2") before shadow-evaporation. Ap- 
proximately 40 nm was removed from the top resist layer 
during the ion mill. The junctions were oxidized for 
80 mins at 5.0 mBar. Lift-off was done in N-methyl- 
2-pyrrolidinone at 80 °C. 

We find that the junctions age very little, the resistance 
value changes less than 1% over a period of ten days. An 
optical micrograph of the devices is shown in Fig. [STl 
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FIG. SI: (Color online) Optical micrograph of the six 'Xmon' 
qubits on two chips, formed by the Al superconducting (light) 
and the exposed sapphire substrate (dark). The qubits are 
capacitively coupled to readout resonators, which couple to a 
readout line in a frequency multiplexed readout scheme ^]. 
The central linewidth S and gap width W are varied from 
8,4 /xm to 24,24 /xm. The top three Xmons have a single arm 
length of L = 130 /xm, the bottom three have L — 165 /xm. 
Test resonators provide an independent measurement of the 
quality factor. 



FIG. S2: (Color online) Z pulse shape measurement. The 
qubit is excited with a 7r-pulse with 20 ns Gaussian enve- 
lope, while a rectangular wave is applied on the Z line. The 
pulse sequence is shown in the left inset. We use Zamp = 0.2 
(arb. units) [3]. The right inset is a cross section at AZ = 
(squares) and AZ = Zamp (circles), as indicated by the ar- 
rows (right). The solid lines are fits to partial qubit rotations 
from a 7r-pulse with a truncated Gaussian envelope. 



Z PULSE 



In order to quantify the response time of the qubit 
to a Z pulse, we simultaneously apply a 7r-pulse with a 
Gaussian envelope on the XY control and a rectangular 
wave pattern on the Z control line. The pulse sequence 
is shown in the left inset in Fig. [S2l We slide the rectan- 
gular wave in time (Ar) and offset (AZ), while retaining 
the amplitude of the wave Zamp constant. The excited 



state population is plotted in Fig. 
time and amplitude. 



as a function of 



A cross section of the main figure at AZ = (squares) 
and AZ = Zamp (circles) is shown in the right inset. The 
measured response can be accurately described by qubit 
rotation from a partial 7r-pulse (yr-pulse duration: 20 ns), 
with a truncated Gaussian as envelope (solid lines). No 
other time constants are included. We find that the rise 
at AZ = Zamp is best described when assuming a 0.5 ns 
delay compared to the fall at AZ = 0. We conclude that 
the qubit frequency is tuned to the desired frequency on 
a timescale of nanoseconds. 

We find the Z control cross-talk between adjacent 
Xmons to be 1.0-1.5 %. 



QUBIT DECAY RATE 



The frequency-dependent decay rate for the six qubits, 
shown in Fig. [ST| is displayed in Fig. [S3l 
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FIG. S3: Frequency dependence of Fi for six qubits with 
different S and W. The frequency stepsize is 5 MHz for S,W 
= 8,4 /xm and 2 MHz otherwise. 
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FIG. S4: (Color online) Qubit excited state probability versus 
time, in a simulation of a qubit-defect system. The qubit is 
on resonance with the defect, whose energy relaxation time is 
Fi^D- Here, we used a coupling strength g/27T — 1 MHz. 



SIMULATION OF A QUBIT-DEFECT SYSTEM 

In order to elucidate the crossover from coherent to 
incoherent decay of the qubit state, we numerically sim- 
ulate a qubit-defect system. The system consists of two 
coupled two- level systems, with coupling strength g and 
defect energy decay rate Fi^d; the qubit is placed on res- 
onance with the defect. The qubit excited state probabil- 
ity is shown in Fig. [S4l For Fi^d < 9^ the excitation co- 
herently swaps back and forth between qubit and defect, 
decaying slowly. When the decay rate exceeds the cou- 
pling strength (Fi^d > g) coherent swapping vanishes and 
an exponential decay appears, as the qubit state decays 
incoherently. Interestingly, the excitation decays most 
quickly for Fi^d ^ ^• 



ANALYTICAL EXPRESSION FOR LOSS IN A 
QUBIT-DEFECT SYSTEM 

Here we derive an analytical expression for the energy 
loss rate arising from a qubit coupling to a single two-level 
defect. We consider a system with two coupled two- level 
systems. We solve the master equation in the Lindblad 
form, 



^[i7,p] + ^D[a]p, 



(SI) 



with V[C]p = CpC+ - (C+Cp + pC+C)/2 and the 
Hamiltonian given by 

H = fyjjQO^a + huj^b'^b + %(6+ (g) a + 6 (g) a+), (S2) 

with a and b the lowering operator for qubit 
and defect, respectively, and u the radial tran- 
sition frequency. We model Markovian decoher- 
ence through the Lindblad terms: Ci=i_4 = 
{ay/Fi^Q, a^a^/2^^, b^/[\^, b^b^/2^^}, denoting 



energy and phase relaxation for qubit and defect, respec- 
tively. Here we used the number operators a'^a and b'^b 
to express pure dephasing, and Fi q denotes qubit relax- 
ation. We take Fi^q <^ Fi^d- As we are interested in 
the relaxation of a single excitation, we only consider the 
states {|00), |01), |10)}. In the interaction picture and 
matrix form the above becomes 
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with p the density matrix, and A = ujt — ^q- 

We are interested in the decay of the qubit excited 
state probability /)22, the relevant equations extracted 
from above are 



P22 = 


-^g[P23 - P23) - ri,QP22 


(S5a) 


P23 = 


-W{P33 - P22) + ^Ap23 - rp23 


(S5b) 


Pss = 


-ig{p23 - P23) - Fi,dP33 


(S5c) 



with r = ri,D/2 + r^,D + ri,Q/2 + r^,Q. 

In the limit Fi^d > 5, /O33 ~ 0, and we can approximate 
the system with two coupled differential equations. In- 
serting P22 = exp(-rit) and p2z = {Pr + iPi) exp(-rit) 
gives 



Ti = 2gPi + Ti.Q 



(S6a) 
(S6b) 
(S6c) 



The solution for the qubit energy decay rate Fi in the 
presence of a two-level defect is (for Fi^d > 5 > Ti^q) 



Ti 



2g'r 

F2+A2 



ri,Q. 



(S7) 



MONTE CARLO SIMULATION OF DEFECTS IN 
THE XMON QUBIT 

In order to quantitatively understand the Xmon 
qubit 's energy decay as well as its variation over fre- 
quency, we have performed a Monte Carlo simulation 
for defects in the capacitor. The defect density for 
AlOa^ in tunnel barriers has been established in mea- 
surements with phase qubits [4], with the distribution 
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FIG. S5: (Color online) Monte Carlo simulation for defects 
in the Xmon capacitor, with 5'=8, VK=8 /im and the Xmon 
arm length L=165 /xm. (a) Distribution of the defect coupling 
strength. A total of 210 defects/GHz have a coupling strength 
g/27v > 0.1 MHz. A total of 30 defects/GHz have a coupling 
strength g/27i > 0.22 MHz. (b) Distribution of the defects 
along position; number of defects in each 100 nm wide section 
is plotted. The cross section of the capacitor is dashed, (c) 
The simulated decay rate in a 0.5 GHz bandwidth using the 
distribution in (a,b), and Eg.lSTl 



over dipole moment p, volume, and frequency given by 
po\/l — P^/Pmax/P? with po ^ 10^//im^/GHz, and the 
maximum dipole moment Pmax = 6 D. As the capaci- 
tor metal oxide and exposed substrate both consist of Al 
oxide, we assume that these numbers are a fair represen- 
tation of the defect density in the qubit capacitor. 

We randomly place defects, with the proper distribu- 
tion over p, in a 3 nm thick dielectric layer (e^ = 10) 
on the substrate-air and metal-air interfaces (top metal 
surface as well as the etched edges). The substrate- metal 
interface is assumed to be thoroughly cleaned [1] and to 
contain no significant defect density. Using a thickness 
of 2 nm instead of 3 nm for the dielectric layer does not 
significantly infiuence the results. The coupling strength 
g = pE is then calculated using a simulation for the 
electric fields in our geometry. The results are shown in 

Fig. EH 

In Fig. [S5k . we have plotted the defect distribution 
over coupling strength. We find that, in a 1 GHz band, 
30 defects have a coupling strength of g/27T ^0.2 MHz 
(square). This simulated value is close to the experimen- 
tally observed density of '^30/GHz. These strongly cou- 
pled defects are predominantly located within a '^100 nm 
distance from the etched metal edges, including the ex- 
posed substrate surface close to the metal edges and ca- 
pacitor metal oxide, where the electric fields are largest, 
see Fig. [S5b . 

The simulated qubit decay rate for the same defect dis- 
tribution is shown in Fig. [S5b in a 0.5 GHz bandwidth. 
For the defect decay rate we have assumed a uniform dis- 
tribution based on the experimentally determined values: 
1/F2,D = 50 - 100 ns (F2,D = Fi,d/2 + F^^d). The sim- 



ulated decay rate reproduces both the peaks as well as 
the variation in the background which are observed in 
the measurement. 
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